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We study the quantum Hall transition using the density-density correlation function. We show 
that in the limit h the electron density moves along the percolating trajectories, undergoing 
normal diffusion. The localization exponent coincides with its percolation value u = ^. The frame- 
work provides a natural way to study the renormalization group flow from percolation to quantum 
Hall transition. We also confirm numerically that the critical conductivity of a classical limit of 
quantum Hall transition is a^x = \/3/4. 

I. INTRODUCTION 

Integer quantum Hall transition is a localization-delocalization transition experienced by a quantum mechanical 
particle moving in a two dimensional plane in a magnetic field perpendicular to the plane and in a random scalar 
potential V{x). As the energy of the particle approaches its critical value, the particle becomes more and more 
delocalized. In spite of the fact that the problem has been studied numerically and experimentally for 20 years, it has 
resisted exact theoretical descriptions. 

In this paper we examine a new method to probe the properties of the quantum Hall transitions recently introduced 
by Girvin ai We show that this method provides a natural framework in which the classical limit of quantum 
Hall transition (classical percolation) can be studied. We show, both numerically and analytically, that even in this 
classical limit the percolating particle exhibits normal diffusion, just as the quantum critical particle would. The 
critical exponent ly turns out to be equal to the percolation exponent | . In this framework, the renormalization group 
flow from percolation to the quantum Hall transition can be studied. 

The method introduced in is based on the introduction of the density-density correlation function 

G(x,t) =^(n|p(x,Op(0)|n) (1) 

n 

where p is the density operator and the sum goes over the Hilbert space of states. This correlation function has a 
simple meaning. It is the probability that a quantum particle, localized at the origin at the initial moment in time, 
is found at the position a; at a time t. Indeed, set up a wave function of a particle localized at the origin, 

i;ix) ^Y.^n{0)M^), (2) 

n 

where ijjn are the eigenstates of the Hamiltonian of the problem. At time t the wave function will be (we set ^ = 1) 

^(a:,i)=5]C(0)V'„(a:)e^^"*, (3) 

n 

where En are the eigenenergies. The probability is equal to 
The matrix elements of the density operator are 

Pm,n{x) = Vm(a;)V'«(a;)- (5) 

It is therefore obvious that the correlation function 

G(x,i)=E 

{x,t)pn,m(0) (6) 
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is equal to the probability 

It sounds reasonable that such a correlation function, which tests how far a particle can go in a given amount 
of time, is useful when studying localization. In the localization literature, however, a different function is usually 
studied. It is called the spectral function and is given by the square of the single particle Green's function ||] 

with e and 5 infinitesimally small numbers. As is well known, the conductance of the system is simply related to the 
spectral function S{x,ui; E). This function contains more information than the correlation G. It tells us how far a 
particle can go if it propagates along the state with a given energy E. The correlation function G is related to the 
integral of S over E, 

/oo 
dE S{x,u;E), (8) 
-oo 

where G{x,lu) is a Fourier transform of G{x,t). The last formula can be obtained by directly integrating (jj) and (0). 

One could argue that the density-density correlator G{x, t) we just introduced would not give much information 
about the localization properties. Indeed, in a transition such as Anderson transition in 3D, this function will always 
be long ranged due to the contributions of the delocalized states to (||). In contrast, the spectral function S{x,u};E) 
tests only the states at the energy E. 

The insight of Girvin at al was that in quantum Hall transition only the centers of the Landau bands are believed 
to be delocalized. Moreover, the physics of quantum Hall transitions allows us to restrict the particle to the lowest 
Landau level. In practical terms, that means restrictingthe sum in (|^) to go only over the states in the lowest Landau 
level. Then there is only one delocalized state. In (||) only the states whose localization length is larger than x 
contribute to the integral over all the energies. Consequently the density-density correlation function includes in itself 
some useful information about quantum Hall transitions. 

It follows from the scaling theory of phase transitions that the localization length £,{E) is divergent as we approach 
the localization transition with the exponent 

aE)^\E-E,\-' (9) 

where Ec is the critical energy. It has also been argued Q that a quantum particle in 2D at the localization- 
delocalization transition undergoes normal diffusion, which in particular results in a finite critical conductivity. Heuris- 
tically people sometimes estimate the Fourier-transformed density-density correlator by an interpolating method based 
on (||). Until the particle hits the localization length, we assume that 

S{k, t; E) cx exp (~Dt^^ , if < S,'^{E), (10) 

which describes diffusion with the diffusion constant D. This assumption is true only if A: ^ [Q. The diffusion 
at small k reflects the existence of critical conductivity in quantum Hall effect. On the other hand, at larger k the 
motion of the particle becomes more complicated. That reflects the fact that critical conductance is not Ohmic and 
in fact dependent on the shape of the sample. Obviously, (^0|) will not work if the particle is deeply in the localized 
regime, when \E — Ec\ '> £'max, with i^max describing the scale where (^ and (^ break down. 
If we wait long enough, 

S{k,t-E)^e^^{-^e{E)^, ifDt:$>eiE). (11) 
Roughly speaking, G(x, t) can be estimated as 

G{k,t) (X dE exp I ^ I , (12) 




where the function being integrated interpolates between the limiting behavior (|T^ and ([T^) . Using (^ and changing 
the variables E = Ec + zt~^ , we can estimate this integral to go as at 
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Gik,t)^t-^J^ "-^^d^ exp (- ,(^5+,..) ) - (13) 

oc 1 - t-^g{k'^t), 

where g is some complicated function. We used that the integral on the second line can be extended to infinity. G{x,t) 
is diffusive at small k (a function of ^) with a prefactor which is t dependent. Knowing the prefactor allows us to 
determine the critical exponent i/. This formula has been tested numerically in ^ and was found to give values of v 
for the quantum Hall transition around the value accepted by the community. 

Imprecise as the derivation of this formula may seem, we will show later in this paper how it can be given a precise 
meaning in the classical limit. 



dz 



1 — exp 



kH 



II. EQUATION OF MOTION 



Now we would like to demonstrate that, as was first shown in [0, the density-density correlation function G(x,t), 
when projected to the lowest Landau level and in the presence of a scalar potential V{x), satisfies a simple equation. 
Choose the basis of states in the lowest Landau level as in (see H and references therein) 



'^n{z, Z*) 



: exp 



ZZ 

w 



(14) 



where, as always, z = x + iy, complex coordinates on the 2D plane and I is the so called magnetic length. Consider 
the Fourier transformed density operator 



Pk = exp 



{kz* + k*z) 



(15) 



where similarly to the complex coordinates k ~ k^ + iky. Let us project this operator to the lowest Landau level, 
following the methods developed in Q. The operator consists of the product of two factors. One of them, exp [f 
is analytic in z and so acts entirely within the lowest Landau level. As such, it does not have to be projected. Indeed, 
if we multiply one of the wave functions ( [l^ ) of the lowest Landau level by this factor, the result can be reexpressed 
as a linear combination of other functions of (|l^. The other factor involves z* and therefore does not act within the 
lowest Landau level. It has to be projected. We note, however, that the operator 



^=''d-z''r 



acts within the lowest Landau level, 



Furthermore, 

d^z (z*)"(z*)'z'"exp 
In other words, 



Zipniz, Z*) = V2Pn Ipn-l 



ZZ 

2P 



d'^z {z*)" exp 



4/2 



2f 



dz 



z exp 



zz 
4P 



( n\ {z*y |to ) = ( n 



d 



1 



^'d-z^-2\ 

This allows us to write down the projected density operator as in 

Pk = exp 









—kz 
2 


exp 


— k z 
2 



(16) 



(17) 



(18) 



(19) 



(20) 



Normally the density operators commute with each other, but the projected density operators do not. Multiplying 
two density operators together and using (Rq) we obtain 
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PkPq = exp 



{k ■ q + ik X q) 



Pk+q, 



(21) 



where k ■ q is the dot and k x q is the cross products of two vectors {k^, ky) and {qx, qy)- 

Finally we note that the Hamiltonian of a particle moving in a lowest Landau level in the presence of the scalar 
potential V{x) is just that potential itself projected to the lowest Landau level, 



(22) 



H= dx V{x)p{x) = / dk V{k)p-k 



Therefore, the equation of motion for the operator p is 



d f 

i^^Pfc = [H,Pk] = / d^q 2isi 



sm I —k X q j V{k — q) exp 



{P-k- q) 



where we restored the h factor. Clearly G(k,t) = tr{pk{t)p-k{0)) satisfies 



ih-^G{k, = J 2* sin (^—k x qj V{k - q) exp 



--^(k'-k-q) 



G{q,t), 



with the boundary condition 



G(fc,0) = tr {p{k)p{-k)) = iVcxp ( --k' 



(23) 



(24) 



(25) 



N being the number of states in the lowest Landau level. Solving this equation allows us in principle to determine 
G{k,t). Therefore, it encodes all the properties of quantum Hall transitions. Unfortunately, while this equation is 
simple to write down, solving it is not a trivial task. 



III. CLASSICAL LIMIT 



Let us consider the classical limit of the quantum Hall transition. The classical limit fi, ^ is equivalent here to 
the limit I since P oc H. In this limit, the equation ( p^ ) becomes 



dG{x,t) _ f 
dt ~ h 



eijdiVdjG{x,t). 



Since P (x h, the planck constant completely drops out. In what follows, we choose the units where h/P — 1. 
The formal solution to this equation is straightforward to obtain: 

G{x,t) = 5{x{t) - x), 

where x{t) is a trajectory satisfying the following equation 



dxjjt) 
dt 



e,,d,V{x{t)). 



(26) 



(27) 



(28) 



Studying the solutions to the equation ( |26| ) reduces to studying the trajectories satisfying (|28 

The tangent to the trajectory x{t) is always perpendicular to the gradient of the potential V{x{t)). Therefore, 
V{x{t)) is conserved along the trajectory. In other words, the particle moves along the equipotential lines, as shown 
on Fig. ||. The value of the potential along the trajectory becomes an important parameter Vq = V{x{t)) which 
can be assigned to each trajectory satisfying (28). In a way, the potential V{x) can be thought of as a mountainous 
landscape with peaks and valleys. As water fills the valleys of the landscape to the height Vq, the equipotential lines 
of Fig |l| are the shapes of the lakes formed in this way. 
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FIG. 1. Equipotential lines of a typical random potential V{x). These are trajectory lines of the solutions to (|28[). The 
potential in this example is generated as a sum of Gaussian shapes located at the nodes of a square lattice with random 
amplitudes. 

In this form the problem we study becomes equivalent to classical percolation Indeed, the potential V{x) is 
random and independent in different points in space. We can consider the value of the potential at a point a; as a 
random variable. Take a particular trajectory x{t) with a particular value of V{x{t)) = Vq. Let us say that the point x 
is "on" if V{x) > Vq and "off' otherwise. Then the trajectory surrounds the clusters of the "on" sites. The trajectory 
becomes the so-called "hull" of the percolation cluster, and the area it surrounds is the percolation cluster itself (or 
the islands in the landscape analogy, for Vq > 0). In the percolation problem, the site is "on" with a probability p. 
In our problem it can be expressed as in 



POO 

p= P{V) dV (29) 

JVn 



'Vo 

where P{V) is the probability for the values of the potential at different points in space (sites). In this way, we 
recover the classical percolation picture of the quantum Hall transitions ||] . The electron density percolates along the 
equipotential lines of the scalar potential. The trajectories are closed. However, as Vq oc (p — ^) approaches zero, the 
trajectories become more and more extended in space. Denoting by ^ the root mean square distance between points 
on the trajectory, we say that the percolation hulls exhibit the delocalization transition, with ^ cx 1q 
Define the area of a cluster as 

A= [ d^x 9{V{x) ~ Vo), (30) 



(x) = 1 when a; > and 9{x) ~ when a; < 0. 

It turns out that the time it takes to go along the boundary of a cluster is just dA/dVo, or 



d'^x S{Vo - V{x)). (31) 



This can be derived directly from ([281). Alternatively, one can argue that ( |28D are the Hamiltonian equations, with V 
being the Hamiltonian. Then S is the "adiabatic invariant" and its derivative is obviously time. 

From here one can conclude that, since the area of a cluster goes as A cx Vq"'^ , the time goes as T cx Vq ^ ^, where 
7 = yI is another well known percolation exponent However, as we will see below, T is not equal to the time it 
takes for a particle to go around the cluster. Rather it is much bigger having to do with the "swiss cheese" shape of 
a typical cluster. As a result, T sums the time it takes go around the cluster, and around each hole in the cluster. 

The classical percolation picture of the quantum Hall transitions was suggested a long time ago. It was subsequently 
demonstrated that it does not capture the physics of quantum Hall transitions. Indeed, the exponent i' computed 
in this way coincides with the its value for percolation = | 0, while the true value for v is believed to be around 
h' — 2.35 w |. We can say that the classical density-density correlation function G{x) forgets about the phase a 
quantum particle carries with itself. It was also known that the classical percolation picture can be derived as a 
classical limit of the Chalker-Coddington model which is believed to lie in the same universality class as quantum 
Hall transitions. 

Some time ago, F. Evers undertook a thorough numerical study of the solutions to the equations ( p8| ) at zero energy 
Vb, motived by classical limit of quantum Hall transition, in the paper H. Its results will be useful for us below. 
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Now we are going to study percolation as a classical limit of quantum Hall transition. The approach we are using 
allows us to give a precise definition of what percolation is for quantum Hall transition, which quantum Hall quantities 
percolate in the classical picture and finally how they change as we allow for finite I to account for the quantum effect. 



IV. PERCOLATION 

To study the percolation limit, we suggest a lattice model of the equation (p^). As a possible lattice model, we 
could have taken the classical limit Chalker-Coddington model This is the version of the model where a particle 
turns left or right at each node with a certain probability, instead of a certain quantum probability amplitude as in 
the original work However, to make a closer connection with percolation, we would like to consider a somewhat 
different model. 

Consider a triangular lattice. Associate with each of its sites a random number V drawn from a probability 
distribution P{V). Now associate the points of the dual lattice with the path x{t). Each path is assigned a number 
Vq, the value of the energy on this particular path. At each step the path moves in such a way that the potential of 
the left node of the bond the particle crosses is greater than Vq which in turn is greater than the potential of the right 
node. For example, the particle on Fig || enters the triangle through the horizontal bond, as shown by the solid line. 
Therefore, a > Vq > b. If c > Vq, it would leave the triangle through the right bond, as shown by the dashed line. If 
c < Vq, then it would leave the triangle through the left bond. 



c 




FIG. 2. The particle moves in such a way that the potential c of the left node of the bond the particle crosses is greater than 
Vo which in turn is greater than the potential b of the right node. 

The trajectory then becomes the hull of percolation clusters. A typical hull is shown on Fig 0. 




FIG. 3. A typical trajectory. The value of its potential Vo = 0.4 

The equation ( p6[ ) also contains dynamical information. The velocity of the particle which follows (^) is proportional 
to the gradient of the potential. We model it by saying that the time it takes for a particle to traverse a bond of the 
triangular lattice is inversely proportional to the potential difference at the ends of the bond. For example, the time 
it takes to traverse the bond of the Fig. |^ (the process denoted by the dashed line) is l/(c — 6). 

Using this model, we performed extensive numerical simulations. In the simulations we assumed that the values of 
the potential at the nodes of the triangular lattice are random numbers uniformly distributed on the interval [—1,1]. 
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The lattice size was 8500 by 8500, and we discarded the trajectories which went outside the lattice. The typical 
trajectories were, however, of much smaller size, about 500. We discovered that the time it takes for a particle to go 
along a trajectory is essentially proportional to the trajectory's length. This is because the velocity of the trajectory 
remains roughly the same as Vq is varied, having no singularity or zero as Vq passes through 0. Indeed, it is easy to 
estimate the average velocity by assuming that all the potentials at the corners of a triangle are random independent 
numbers. Then all we need to do is to average l/{Vi — V2) where Vi and V2 are the values of the potential at the 
ends of the bond the trajectory is traversing. For the potential uniformly distributed on the interval [—1,1] we obtain 
that the average velocity is very weakly Vq dependent, with its value being log(4) when Vq = 0. This was confirmed 
numerically. In fact, that coincides with the earlier observation in 

Interestingly, that allows us to find the time T it takes for a particle to go around the cluster. The length of a 
trajectory is known to go as 



L cx V„-'^ 



(32) 



where 7' = 2 ||l^. Since the average velocity is only weakly Vq dependent, we obtain that T cx V^^. That should be 
contrasted with the time it takes to go around the cluster and around each hole in the cluster which we found in the 
previous section. 

Replacing time with length we studied x'^{t) where the average means averaging over different trajectories with the 
same value of Vq. In our simulations we averaged over about 10^ trajectories. Fig. ^ summarizes the results of the 
numerics. 



7, , 8000 
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FIG. 4. x'^{t) computed for trajectories with various Vq. The highest line represents the lowest Vq. The potential at the 
nodes of the lattice is uniformly distributed from —1 to 1. Time t is simply the number of steps, is measured in units where 
the triangular lattice size of Fig. ^ is one. 



It is clear that x"^ (t) is a straight line 



x^{t) ^ Dt, i < t* 



until a certain time t* which varies with Vo- After that. 



(33) 



(34) 



with the parameter ^ depending on Vq but not on time. 

The length ^ is none other than the localization (correlation) length of the percolation problem we introduced 
above. Indeed, if we wait long enough, the particle satisfying (Gq) starts winding around the closed trajectories such 
as on Fig. |. Then x'^{t) measures the mean square distance between two points on the percolation cluster's hull 
which is the trajectory we study. So we can write down 



£.^Vq 

which we also confirmed by our numerical simulations. 



(35) 
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The numerical result that x'^{t) oc t ior t < t* shows that the particle undergoes a normal diffusion. While this is 
what is expected of a quantum particle in a quantum Hall transition [^j , it is far less obvious that a classical particle 
should do the same. It is also not obvious how to explain this from a purely percolation point of view. 

This result was previously obtained in where for the trajectories at Vq = was studied. There it was 

shown that the diffusion was due to a "mysterious" cancellation of exponents in percolation. Perhaps, a certain scaling 
relation between those exists which leads to the cancellation, but at the moment we are not aware of it. 

Now we would like to note that a convenient way to compute would be just summing all the trajectories of Fig. 
^ together. As above, we expect this to be 

/ dVo P{Vo) , (36) 

Vo •' Dt + ^0 

Numerically this exponent is around 0.62 which corresponds to ly = 1.30, very close to the exact value. This method 
proves to be very convenient when studying quantum Hall transition in the quantum regime as well . 
Knowing x{t) it is not difficult to compute G{x,t) according to (|27|). 

G{k, t)=Y^ expiikxit)), (37) 

Vo 

SO that the average is taken over all the trajectories. Compare this formula with (p|). We are led to conjecture that 



S{x, t; Vo) = 2ttP{Vo) exp {ikx{t)), (38) 

where the average, as always, is taken over the trajectories with the same Vo- Indeed, not only is automatically 
satisfied, but also this form of S reflects the fact that the particle propagates with a given energy. 

With this definition of S and G it is easy to check that the formulae (pX|)-(p^ are indeed true and precisely defined. 
Indeed, their applicability was at small k where 

fc2 / p \ 

exp{ikx{t)) « 1 - —x^it) w exp i-—x^{t)\ . (39) 

From this point on, (^^-(|l3|) can be carried on without any further obstacles. The only difference between classical 
and quantum case will be the value of ly, which in the classical case is equal to 4/3. 

It is interesting to note that at large k the functions S and G will be very different in quantum and classical case. In 
the quantum case, 5 is a function of the ratio g^/w. In the classical case this question was investigated in [^j. There 
it was shown that S{x, t; 0), computed at zero energy, was a function of /uj, with 7/4 being the fractal dimension of 
the infinite percolating cluster. We refer the reader to |^ for further information on this subject. In the meanwhile, 
the form of G or even S at generic values of energy in the classical case remains unknown. 

Finally, we would like to comment on the interesting feature of the numerical data of Fig. ^. The diffusion constant, 
when read off this data (given in units where the bonds of the triangular lattice of Fig. have unit length), is equal 
to 

D = 1.008 ±0.014 (40) 

This number is remarkably close to 1 

Assuming that _D = 1 is indeed an exact relation, we can demonstrate that the critical conductivity is equal to 



(41) 



4 /i ' 

an exact result recently derived in [p^ by a very different approach. To do that, we need to use the relationship 
between the diffusion constant and the critical conductivity proposed by F. Evers and W. Brenig in |l^ ]. 
According to the Einstein's relation, the conductivity is given by 

(Ta^x = e^Pj, (42) 

where e is the electron charge and p is the density of states per unit area per unit energy at the transition. The 
problem is to find p in the lattice model introduced here. 
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We can use a standard quasiclassical argument to find the density of states for a particle moving in a magnetic field 
and a potential V{x). The total number of states whose energy is below E is given, in the quasiclassical limit, by the 
area of the region where V{x) < E, divided by 2-^1"^, 



A{E) = ^ j d'xe{E-V{x)). (43) 



The density of states per unit area is 

2ttPA 

where A is the total area. 



p{E) = -^ I d'xS{E~V{x)), (44) 



Comparing p{E) with (31) we see that p{E) coincides with the time T it takes for a particle to go over all the 
trajectories with energy E. Restoring h and I dependence in we find 

ME) ^ ^ (45) 

Therefore, the conductivity reduces to 

It is now clear that the overall answer is independent of the units in which we measure D, as long as they are the 
same as the ones for A and T. 

In the lattice model of Fig. ^ the total time is equal to the total number of links available for a particle with zero 
energy, the energy of the transition. Assume that the total number of triangles is equal to N. The area A of the 
lattice is equal to iV x where ^ is the area of a single triangle. The total number of links is L = We took 
into account that each triangle has 2 links with a probability 3/4. This is because with a probability of 1/4 all three 
energies at the vertices of a triangle are either less or greater than 0, in which case there are no trajectories in the 
triangle at all. But each link belongs to two triangles, hence the overcounting and the answer of Substituting 



this into (46) and taking into account that « 1, we arrive at 



<^xx — L) — « — , (47) 

in very good agreement with jl^ . 

This result is also in agreement with the one of where the critical conductivity was estimated at (Jxx ~ 0.45^. 

Notice that the conductance of a finite sample would not be Ohmic since the particle does not undergo a simple 
random walk, as reflected in a complicated structure of G or S' at finite momenta k. Only for very thin samples, 
when it takes little time for a particle to reach one end of a sample from another, can the Ohm's law be applicable, 
in agreement with |p^ . 



V. CONCLUSIONS 



To summarize, we showed that in the classical limit the density of a particle undergoing the quantum Hall transition 
moves along the percolation hull trajectories. At small times this particle exhibits normal diffusion, while at large times 
it localizes with the localization length exponent given by the percolation v = ^. Our paper sets a new framework in 
which the classical percolation to quantum Hall effect renormalization group flow can be studied. Finally, we confirm 
numerically with a high degree of accuracy the result a^x = \/3/4 for the critical conductivity of a classical quantum 
Hall transition. 
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